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ABSTRACT 


This  report  examines  the  sensitivity  of  material  response 
calculations  to  the  choice  of  equation  of  state  model.  Three 
equation  of  state  models,  all  of  which  are  available  as  sub¬ 
routines  in  the  RIP  code,  are  considered:  1)  modified  PUFF; 

2)  RIP  mixed  phase;  and  3)  GRAY.  Each  of  these  models  is 
used  to  calculate  the  one-dimensional  response  of  aluminum, 
beryllium  and  titanium;  this  series  of  calculations  is  per¬ 
formed  for  two  distinct  x-ray  sources.  The  sensitivity  of 
the  calculated  material  response  to  the  choice  of  equation  of 
state  model  is  characterized  in  terms  of  the  generated  impulse 
and  the  peak  propagating  stress  at  the  time  the  radiation  source 
is  cut  off.  For  the  calculations  presented  in  this  report,  the 
three  equation  of  state  models  are  in  fairly  good  agreement. 
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1.  INTRODUCTION 


1.1  OBJECTIVE  AND  PROCEDURE 


The  purpose  of  this  report  is  to  determine  the  sensi¬ 
tivity  of  material  response  calculations  to  the  choice  of 
equation  of  state  model.  Three  equation  of  state  ..models ,  all 
of  which  are  available  as  subroutines  in  the  RIP1  J  code,  are 
considered:  1)  modified  PUFF;  2)  RIP  mixed  phase;  3)  GRAY. 

Each  of  these  models  is  used  to  calculate  the  one-dimensional 
response  of  aluminum,  beryllium  and  titanium;  this  series  of 
calculations  is  performed  for  two  distinct  x-ray  sources.  The 
sensitivity  of  the  calculated  material  response  to  the  choice 
of  equation  of  state  model  is  characterized  in  terms  of  the 
generated  impulse  and  the  peak  propagating  stress  at  the  time 
the  radiation  source  is  cut  off.  Whenever  possible,  handbook 
values  are  used  for  the  required  equation  of  state  data;  no 
attempt  was  made  to  adjust  the  data  in  order  to  reduce  the 
differences  between  the  three  models. 


1.2  EQUATION  OF  STATE  MODELS 

The  first  equation  of  state  considered  in  this  report 
is  referred  to  as  the  modified  PUFF  equation  of  state.  I ^ 

The  modified  PUFF  equation  of  state  is  the  simplest  of  the 
three  models;  it  consists  of  a  Mie-Gruneisen  solid  equation 
of  state  for  the  compressed  region  which  is  joined  at  the 
ambient  density  to  a  function  which  limits  to  an  ideal  gas 
equation  of  state  as  the  density  approaches  zero.  Pressure 
continuity  is  maintained  at  the  ambient  density  for  all  ener¬ 
gies,  but  the  sound  speed  is  continuous  only  at  the  zero 
reference  energy;  the  discontinuity  in  sound  speed  becomes 
more  severe  as  the  energy  is  increased. 

Ths  RIP  mixed  phase  equation  of  state  is  based  on  a 
virial  expansion  which  was  applied  by  Tuerpe,  Keeler  and 
McCarthy  l  *1  to  the  study  of  aluminum  at  high  temperature  and 
pressure.  The  equation  of  state  was  later  generalized  to 
other  metals  by  D.  Steinberg  and  reported  in  an  internal  memo 
dated  August  22,  1967  at  Lawrence  Livermore  Laboratory. 
Incorporation  of  the  Steinberg  equation  of  state  into  the  RIP 
code  involved  the  following  tasks: 

1.  The  caloric  equation  of  state  was  inverted 
and  the  resulting  expression  for  the  tem¬ 
perature  was  substituted  into  the  thermal 
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equation  of  state,  thereby  yielding  a  formula 
for  the  pressure  as  a  function  of  energy  and 
density.  Seven  independent  coefficients 
appear  in  this  formula  and  they  are  evaluated 
by  the  EOS  GEN*  code  and  subsequently  input  to 
EOSTAB  common  in  the  RIP  code. 

2.  A  Maxwell  construction  was  required  in  order 
to  determine  corresponding  values  of  pressure, 
energy,  sound  speed  and  density  on  the  mixed 
phase  boundary.  This  task  is  also  performed 
by  the  EOSGEN  code,  and  points  on  the  mixed 
phase  boundary  are  subsequently  input  to  EOSTAB 
common  in  the  RIP  code. 

3.  A  join  function  was  developed  which  provides 
for  pressure  and  sound  speed  continuity  between 
the  virial  expansion  for  expanded  states  and 
the  Mie-Gruneisen  solid  equation  of  state  for 
the  compressed  region.  The  coefficients  in  the 
join  function  are  evaluated  by  the  EOSGEN  code 
and  subsequently  input  to  EQST  common  in  the 
RIP  code. 

A  brief  discussion  of  the  EOSGEN  code  and  associated 
corrections**  to  the  E0S2  subroutine  in  the  RIP  code  is  pre¬ 
sented  in  Appendix  A.  All  of  the  work  presented  in  Appendix  A 
was  performed  in  support  of  Contract  DAAD05-73-C-0489  and  was 
required  in  order  to  perform  hydrodynamic  calculations  using 
the  RIP  mixed  phase  equation  of  state  in  the  RIP  code. 

The  third  equation  of  state  considered  in  this  research 
is  the  GRAY  three-phase  equation  of  state  for  metals  developed 
by  E.B.  Roycef3]  at  Lawrence  Livermore  Laboratory.  Material 
near  normal  density,  in  the  solid-liquid  region,  is  described 
by  a  scaling- law  equation  of  state  for  metals  developed  by 
Grover. I* J  The  scaling-law  equation  of  state  includes  a 
Gruneicen  description  of  the  solid! normalized  to  the  experi 
mental  Hugoniot.  Material  in  the  liquid-vapor  region  is 
treated  according  to  an  equation  of  state  developed  by  Young 
and  Alder;  [^1  their  van  der  Waals  model  uses  an  analytic  repre 
sentation  of  the  classical  hard-sphere  equation  of  state  with 
a  van  der  Waals  attractive  term  added.  The  complete  equation 


Documentation  of  the  EOSGEN  code  is  in  progress. 

It  was  necessary  to  develop  a  new  join  function  in  order  to 
allow  pressure  and  sound  speed  continuity  between  the  com¬ 
pressed  and  expanded  equations  of  state. 
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of  state,  as  described  In  Reference  3,  Is  developed  by  analytically 
joining  the  Grover  scaling  law  and  the  Young-Alder  models  at  a 
volume  In  the  range  1.3  to  1.5  times  normal  volume.  This  Is  accom¬ 
plished  by  adding  correction  terms  to  the  Young-Alder  equation  of 
siate  employed  on  the  low  density  side  of  the  join  volume. 

It  was  necessary  to  make  several  corrections  to  the 
E0S3  subroutine  In  the  RIP  code  before  performing  calculations 
with  the  GRAY  equation  of  state.  These  corrections  are  presented 
in  Appendix  B. 

1.3  SOURCE  DESCRIPTIONS 


The  first  radiation  source  (SPEC  1)  was  prescribed  by 
the  Ballistic  Research  Laboratory. [7D  The  energy  spectrum  and  the 
time  dependence  of  the  energy  flux  are  presented  in  Figure  1  and 
Figure  2  respectively.  The  second  peak  in  the  energy  spectrum  Is 
due  to  a  typographical  error  In  Reference  7.  The  calculations  were 
performed  before  this  error  was  recognized,  but  since  the  same 
spectrum  was  used  throughout,  the  calculations  still  represent  a 
valid  sensitivity  study.  A  fluence  level  of  5952  cal/cm2  was  used 
for  all  calculations  performed  with  this  radiation  source. 

The  second  radiation  source  (SPEC  2)  considered  In  this 
study  Is  a  2-keV  blackbody  with  a  10-nsec  square  pulse  shine 
time.  Appropriate  fluence  levels  were  selected  in  order  to 
guarantee  that  several  zones  of  material  are  vaporized  at  the 
front  surface  of  each  metal  considered.  The  fluence  levels 
selected  are: 


Aluminum: 

Beryllium: 

Titanium: 


100  cal/cm^ 
2000  cal/cm2 
100  cal/cm2 


Normalized  intensity  (10“3  cal/cm2AeV) 


r  I 


Figure  1.  Incident  energy  spectrum  for  SPEC  1  radiation 
source  (normalized  to  unity)  . 
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Normalized  energ; 


Figure  2.  Normalized  energy  flux  vs  time  for  SPEC  1  radia 
tion  source. 


2 . 1  MODIFIED  PUFF 


The  modified  PUFF  equation  of  state  parameters  for 
aluminum,  beryllium,  and  titanium  are  presented  in  Tables  1, 

2,  and  3  respectively.  All  of  the  parameters  appearing  in 
the  tables  were  obtained  from  Reference  8  with  the  exceptions 
of  E lv,  E^,  and  ESM  which  were  calculated  as  described  in 
the  following  paragraph. 

The  incipient  melt  energy,  Eg*,,  is  obtained  by  multi¬ 
plying  the  difference  between  the  melt  temperature  (Refer¬ 
ence  9,  page  D-103)  and  the  reference  temperature,  20°  C,  by 
the  heat  capacity  at  constant  pressure  (Reference  9,  D-97)  . 
The  complete  melt  energy.  Elm,  is  obtained  by  adding  the  heat 
of  fusion  (Reference  10)  to  the  incipient  melt  energy. 

Finally,  an  approximate  value  for  the  incipient  vaporization 
energy,  Elv»  is  obtained  by  multiplying  the  difference  be¬ 
tween  the  melt  temperature  and  the  boiling  temperature  (Refer¬ 
ence  9,  page  D-103)  by  the  heat  capacity  and  adding  the  result 
to  the  complete  melt  energy.  These  three  energy  parameters 
are  used  for  the  momentum  edits  (the  momentum  of  each  thermo¬ 
dynamic  state  is  displayed)  and  for  determining  the  energy 
density,  ESm»  at  which  the  material  should  have  zero  strength. 
All  of  the  data  for  modified  PUFF  is  stored  in  the  ESTCON 
array  in  EQST  common. 


2.2  RIP  MIXED  PHASE 

The  first  thirteen  parameters  listed  in  Tables  1,  2, 
and  3  are  required  as  input  to  EQST  common  in  order  to  apply 
the  EOS2  equation  of  state  subroutine  in  hydrodynamic  calcula¬ 
tions  on  aluminum,  beryllium  and  titanium,  respectively.  The 
remaining  equation  of  state  data  are  generated  by  the  EOSGEN 
code  which  is  described  in  Appendix  A. 

Fundamental  material  properties  of  aluminum,  beryllium, 
and  titanium  which  are  required  as  input  data  to  the  EOSGEN 
code  are  presented  in  Tables  4,  5,  and  6  respectively.  The 
ambient  densities  are  the  same  values  used  for  the  modified 
PUFF  equation  of  state  and  the  reference  temperature  for  all 
three  materials  is  assumed  to  be  300°  K.  The  bulk  modulus 
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of  aluminum  was  calculated  according  to  the  relation 


where  is  the  bulk  modulus,  Y  is  Young's  modulus  and  o  is 
Poisson^s  ratio.  Numerical  values  for  Y  and  o  were  obtained 
from  Table  3.1  on  page  44  of  Reference  1.  The  bulk  moduli 
of  beryllium  and  titanium  were  taken  equal  to  the  linear 
Hugoniot  coefficients,  Ai,  used  in  the  modified  PUFF  equation 
of  state.  Constant  pressure  heat  capacities  obtained  from 
page  D-97  of  Reference  9  were  used  for  cvo»  the  constant 
volume  heat  capacity  at  the  reference  state;  heat  capacities 
at  constant  pressure  and  volume  are  almost  identical  for 
solids  at  low  temperatures.  The  critical  pressure,  tempera¬ 
ture,  and  volume  of  aluminum  were  obtained  from  Table  3.1  on 
page  44  of  Reference  1.  The  critical  points  of  beryllium  and 
titanium  were  calculated  using  a  procedure  presented  in  Refer¬ 
ence  6.  The  required  relations  are 

vc  =  (2.417xl024)  a3  , 

.7232  a 


where  Vc,  Tc,  and  Pc  are  the  critical  volume,  temperature,  and 
pressure,  respectively,  a  is  the  hard  sphere  diameter,  ay  is 
the  coefficient  of  the  attractive  potential  for  the  vapor,  and 
R  is  the  universal  gas  constant.  The  hard  sphere  diameter  is 
related  to  the  excluded  volume,  Vb,  by  the  equation 


where  N  is  Avogadro's  number.  Numerical  values  of  ay  and  Vfc 
are  presented  in  Table  3.2  on  page  52  of  Reference  1. 
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2 . 3  GRAY 


The  GRAY  equation  of  state  parameters  for  aluminum, 
beryllium,  and  titanium  are  presented  in  Tables  11,  12  and 
13  respectively;  the  parameters  also  used  in  the  modified 
PUFF  equation  of  state  are  not  repeated  in  these  tables. 

The  entropies  of  melting  were  obtained  from  Reference  9;  all 
other  GRAY  equation  of  state  parameters  were  obtained  from 
the  original  Lawrence  Livermore  Laboratories  report  by 
E.B.  Royce . 15 j 
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MODIFIED  PUFF  PARAMETERS  FOR  ALUMINUM 


MODIFIED  PUFF  PARAMETERS  FOR  TITANIUM 
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EOSGEN  INPUT  DATA  FOR  ALUMINUM 
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EOSGEN  INPUT  DATA  FOR  TITANIUM 
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TABLE 
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gm/cm  -2.700  3.834  8.517 


TABLE  8 


/EOSTAB/COMMON 


Purpose :  Stores  the  constants  and  the  tabular  values  which 

describe  the  mixed  phase  region  in  the  RIP  mixed- 
phase  equation  of  state  model. 

used  In:  CARDS,  DREAD,  EDIT,  EOS2,  RIP,  SPRINT 


Q|  g 

Location 

Definition 

THETA (I) 
1=1,100 

Y1 

Temperature  (°K)  of  the  points  for 
the  mixed  phase  region  table.  List¬ 
ed  in  ascending  order  from  T  =  0  to 

T  =  T  ...  These  values  are  deter- 
cnt 

mined  by  the  EOSGEN  code  and  are  not 

used  in  subroutine  EOS2. 

PU(I) 

1=1,100 

Y101 

2 

Pressure  (dynes/cm  )  of  the  points  in 
the  mixed  phase  region  table.  These 

values  correspond  to  THETA  values  and 

are  in  ascending  order  from  P  =  0  to 

P  =  P  ... 
cnt 

RHOU(I) 

1=1,100 

Y201 

Density  (g/cm  )  of  the  points  for  the 

range  p_  >  p  >  p„..  of  the  mixed 
co  cnt 

phase  region:  the  variable  p  is 

CO 

the  value  of  the  density  at  the  cut¬ 
off  point  for  the  mixed  phase  region 
which  is  entered  into  ESTCON (16 ,N)  . 

The  variable  Pcr^t  corresponds  to  the 
density  at  pcr^t*  These  values  corre¬ 
spond  to  the  values  appearing  in  the 

PU (I)  array. 

EU  (I) 

1=1,100 

Y301 

Specific  internal  energy  (erg/g) 
corresponding  to  the  points  in  the 
RHOU(I)  array. 
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TABLE  8,  /EOSTmB/COMMON  (Continued) 


Variable 

Location 

Definition 

CU(I) 

1=1,100 

Y401 

Sound  speed  (cm/sec)  corresponding  to 
the  points  in  the  RHOU(I)  array. 

RHOL(U) 

1=1,100 

Y501 

Density  (g/cm^)  of  the  points  for  the 
range  Pcr^^.  >  P  >  0  of  the  mixed 
phase  region.  These  values  correspond 
to  the  values  appearing  in  the  PU(I) 

array. 

EL  (I ) 

1=1,100 

Y601 

Specific  internal  energy  (erg/g) 
corresponding  to  the  points  in  the 
RHOL(I)  array. 

CL  (I) 

1=1,100 

Y701 

Sound  speed  (cm/sec)  corresponding 
to  the  points  in  the  RHOL(I)  array. 

z  (I) 

1=1,15 

Y801 

Z (1)  through  7  (8)  are  defined  in  the 

EOSGEN  code.  These  are  coefficients 

used  in  the  calculation  of  pressures 
and  sound  speeds  for  state  points 
above  and  below  the  mixed  phase 
region.  Z(9)  through  Z(15)  are  coeffi¬ 
cients  calculated  within  subroutine 

E0S2  and  are  used  for  state  points 
above  and  below  the  mixed  phase. 
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TABLE  9 


EOSTAB  INPUT :  ALUMINUM 


ri  = 

3.0000*02* 

1.8590*03* 

2.0  760*03* 

2.2980*03 

2.7920*03* 

2.9692*03* 

3.1360*03* 

3.9333*33 

3.6520*03* 

9.C79C*r3, 

9.2560*03* 

9.5180*03 

9.9620*03* 

5.1893*33* 

5.9363*03* 

5.6230*33 

6.0720*03* 

6.2990*  03* 

6.5160*03* 

6.7380*03 

7.1820*03* 

7.9093*33* 

7.62  60*23* 

7.8930*23 

8.2920*03* 

8.5190*03* 

8.7260*03* 

8.9560*03 

9.9020*03* 

9.6292*03* 

9.89 60  *C  3 , 

1.0059*39 

1.0512*09* 

1 .1900*09 

1.0739*09* 

1.0356*09* 

1.1178*09 

YlOlr 

0.0033  * 

1.5397*33* 

1.11 92*09* 

5.9833*39 

5. 9992*0'. 

1.5005*06 , 

3.2396*06* 

6.9626*06 

1.9609*07* 

3.1078*37* 

9.6368*07* 

6.7821*37 

1.2820*08* 

1.6902*08, 

2.1787408* 

2.7518*Ce 

9.1719*08* 

5.0253*03, 

5.97  96*08* 

7.0355*33 

9.9619*03* 

1.0833*09, 

1.2310*09* 

1.3895*09 

1.7382*09* 

1  .9283*3  3, 

2.1239*03* 

2.3398*09 

2.7920*09* 

3. 0330*  09, 

3.2838*0!.  • 

3.5999*09 

9.0935*09* 
S. 3000*09 

9.3319*09* 

9.67  92*09, 

9.9853*09 

7231  = 

?. 70CC*QC» 

2.1086*00, 

2.0925*00, 

1.9799*00 

1 .3633*03* 

1.3391*03, 

1  .7570*00, 

1.7371 *30 

1.6131*00* 

1.5687*00, 

1.5257*00, 

1.9893*00 

1.9051*00* 

1.3672*30, 

1.3109*00, 

1.2996*30 

1.2259*00, 

1.1920*00, 

1.1593*00, 

1.1272*00 

1.0697*00* 

1.0391*03, 

1.0390*00* 

9.7915--1 

9.1525-01* 

6. 86C9-C1 , 

8.5689-01* 

B. 2771-01 

7.5885-01* 

7.3886-01* 

7.0327-01* 

5.7635-31 

6.0975-01* 

J.7C37-*! 

5.7260-01* 

5.3090-01* 

9 .8003-01 

7301  = 

0.0003  * 

2.2583*10* 

2.6299*10* 

2.9951 *1C 

3. 7980*10* 

9.1273*10, 

9.5079*10. 

9.6895*10 

S  .*536*13'* 

6.0357*13, 

6.91 73*10* 

6.7997*10 

7.5613*10* 

7.9919* 1C, 

8.3216*10, 

8 . 70«J.S*  1C 

9.9587*10* 

9.8371*10, 

1.3215*11* 

1.0599*11 

1.1350*11* 

1.1729*11* 

1.2109*11* 

1.2990*11 

1.3255*11* 

1.3693*11* 

1  .9328*11  * 

1 .9913*11 

, 

1.5211*11* 

1.5616*11* 

1.6C28*11* 

1.6996*11 

1 .7333*11* 
1.9963*11 

1.7303*11* 

1.3315*11* 

1  .8933*11 

7931  = 

5.0000*05* 

9.3552*05* 

9.2729*05* 

9 .1913*05 

9.3331*05* 

3.9555*05* 

3  .37  87*05* 

3.8029*35 

3.6518*05* 

3.5771*05, 

3.5229*05* 

3.9269*05 

3.2799*05* 

3.2059*35* 

3.131 3*05* 

3.0559*35 

2.5200*03, 
3. .30203. 
9.79CP*n3« 
5. 850303. 
6.96CC03. 
3 .0703 ♦ 03  . 
9. 1 8CC*  C3» 
1.0292  03* 


2  .323305* 
1.1665*07, 
9.9657*07  . 
3.9196*08* 
3.196303. 
l. 55664*9, 
2.562909. 
3.8133+rb. 


1.92C2+CC. 
1.659200* 
1.333COC* 
1.259600* 
1.0957400. 
3.3361 -01 » 
7.  9891-01 . 
6.3323-31. 


3.3737410, 
5.2715410. 
7.1^C°4lO» 
9.C8CC41C, 
1 .0972.11 » 
1.2871411, 
1  .3313*11 . 
1.6861411, 


3.1116*05, 
3.725305. 
3.35:,j4',5. 
2  .39210  5, 
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TABLE  9,  EOSTAB  INPUT:  ALUMINUM  (Continued) 


f53l: 


Y6C1  r 


T73i  = 


YBOlr 


2.937  0*05 
2.5249*05 
2.1233*25 
1. 6876*05 
I .1939*05 
7.  8535*04 

o.ooco 

5.8130-35 
1. 3646-03 
7. 1401-03 
2.0069-02 
1.1567-02 
7. 3595-02 
1.1965-21 
1. 9213-01 
3.7037-01 

1.3191*11 

1.5391*11 

1.5378*11 

1.7323*11 

1.3223*11 

1.9001*11 

1.9712*11 

2.0736*11 

2.2651*11 

1.9969*11 

3.3612*04 
1.2396*05 
1.4227*05 
1 .5720*05 
1.6535*0' 
1.5679*05 
1.6140*05 
1.4814*05 
1.2335*05 
7.8535*04 

3.7658*06 

9.0000*06 

-5.4376*10 


2.8316*35 

2.4464*05 

2.0395*35 

1.5942*05 

1.0839*35 


2.1913-07 
1.3459-  34 
2.C534-C3 
9.0930-  33 
2.36E4-C2 
4.7153-32 
8. 1479-02 
1.3129-31 
2. 1336-01 


1 .4592*11 
1.5546*11 
1.6571*11 
1.75C5* 11 
1.8370*11 
1.9153*11 
1.9841*11 
2.C3«6*U 
2.0573*11 


9.9676*04 
1.2573*35 
1.4578*05 
1.5937*05 
1.6617*05 
1  .6628*35 
1 . 5  °4  4*  05 
1.4433*35 
1.1605*05 


-2.0254*35 
1.2924*  11 
-  2.4900*10 


2 .7558*05 

2.3  £72*05 
1 .9542*35 

1.4  9ez*cs 
9.6723*04 


1.4248-06 
2  .7439-34 
2.9516-03 
1.1  3  54-32 
2. 7595-02 
6.3361-32 
8. 9  958-02 
1 .4435-01 
2.3957-01 


1.47  92*1  1 
1.5  789*11 
1 .6763*11 
1.7684*11 
1  .85  33*11 
1.9299*11 

1  .9354*11 
2.0482*11 

2  .3554*11 


1.0531*05 
1 .3320*05 

1.4  904*05 
1 .61 27*05 
1.6  673*35 
1 .65  43*05 

1.5  715*05 
1.4302*05 
1.0740*05 


-  9.2422*10 
4.5000*01 

-  8  .55  31  *10 


2.6794  *35 
2.2871*05 
1 .8673*05 
1.3993*05 
3.5334*34 


6.3441-Cf 
5.0555-34 
4,0860-03 
1.3979-32 
3.1932-02 
3.9422-02 
9.9C91-C2 
1  .53.9-31 
2. 7539- Cl 


1  .4932*1 1 
1.5986*11 
1  .6952*11 
1.7860*11 
1 .8613*1 1 
1.9441*11 
2.0332*11 
2.0557*11 
2.0551*11 


1.1078*05 
1.3446*3  5 
1.5203*05 
1 .6233*35 
1.6702*05 
1.6441*35 
1.5453*05 
1 .3515*35 
9.6593*04 


6 .2947*05 

O.OOCO 

6.3584*13 


2 .6025*35* 
2.2069*05. 
1  .7765  *05  » 
1 .29  2C  *C5  . 


2.1341-05. 
9.5857-34* 
5.  474  3-C3* 
1 .5931 -02. 
3.6589-02  , 
S .6253  -02  , 
1.C69G-C.  , 
1.7434-31. 


1.5192*11. 
1.6183*11. 
1.7139*11. 
1.8033*11. 
1 .33'3  *1  1  * 
1.9579*11, 
2.3134*11. 
2.0619*11, 


1.16Cr*C5. 
1  .384  9*05. 
1.5475*° 5. 
1  .5426  *35, 
1.6704*06, 
1 .5335*05, 
1.5153*05, 
1 .2964*05, 


3.2553*13, 
0.C00C  , 
1.1386*1 0 
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TABLE  10 


EOSTAB  INPUT :  BERYLLIUM 


M  = 

3.0000*02. 

2.0760*03* 

2.2980*03* 

2.5206*03* 

2.7420*03* 

4.0740*03* 

4.2960*03* 

4.5180*03* 

4.74  C0*03 • 

4.9620*03. 

5.1840*03* 

5 .4  06C*0  3* 

5.6280*03* 

E .8  5  CO*  C 3  * 

6  .C  720*03. 

S. 294003* 

6.5180*33* 

6.7380*03* 

6.9530*33* 

7  .1823*03. 

7.9040403* 

7.6260*03* 

7.8480*03. 

8 .0  7CC*  0 j  * 

8.2926*03. 

8.5140*03* 

8.7360*03* 

8.9580*03* 

9.1833*33* 

9.4020*33. 

9.6240*03* 
1.0734*04* 
1. 1400*04 

9.8460*03* 

1.0956*01* 

1.0668*04* 

1.1178*04* 

1.0290*04* 

1.0512*04. 

rm  = 

O.OOCC  • 

3.4925*03* 

1.9672*04* 

7.9035*04* 

2.6176*05. 

7.2807*05* 

1.7635*06* 

3.8206*06* 

7.5545*36* 

1.3853*07. 

2. 3831*07* 

3.8866*07* 

6.0556*07* 

9.0727*07* 

1.3139*Ce. 

1.8475*08* 

2.5312*38* 

3.3399*08* 

4.4493*33* 

5.7363*33. 

7.2773*08* 

9.1011*08* 

1.1235*09* 

1.3  7C6*C9» 

1.6549*09; 

1 .9739*09* 

2.3454*39* 

2.75  74*09* 

3.2133*39* 

3.7333*09* 

4.2973*09* 

4.9212*09* 

5.6660*09* 

6.3540*09* 

7.1684*09* 

8.0526*09* 

1.3624*10* 

1 • 8000*10 

9.0097*09* 

1.4990*10* 

1.0342*10* 

1.6448*10* 

1.1153*10* 

1 .234  5  *1 U* 

Y20I  - 

1.8510*00* 

1.7476*03* 

1.7337*00* 

1.7195*30* 

1 .7051*30* 

1.6903*00* 

1.6753*00* 

1.6En«CG* 

1.6445*00* 

1.6285  *C0» 

1.5123*00* 

1.5957*00* 

1.5737*00* 

1  .5513*30* 

1.5435*30* 

1.5254*00* 

1.5068*00* 

1.4878*00* 

1.4683*00* 

1.4483*00* 

1.4278*00* 

1.4067*00* 

1.3351*00* 

1.3629*30* 

1.3399*00* 

1.3164*00* 

1.2921*00* 

1.2E69*OC» 

1  .?410*0C* 

1.2141*00, 

1.1883*00* 

1.1572*03* 

1.1270*00* 

1.0954*30* 

1  .0622*00. 

T331  = 

1.0272*00* 

8.3834-01* 

4.4261-01 

9.9011-01* 
7 .4667*01  * 

9.5  644- 01  • 
6.G873-P1* 

9.0764-01, 

8.6069-01. 

O.OOCC  * 

3.3674*10* 

3.8654*10* 

4  .24  78*  10* 

4.6948*10. 

5.1459*10* 

5.6042*10* 

6.0569*10* 

6.5355*1  0* 

7 .3102*1 n. 

7.4915*10* 

7.9797*10* 

8.4752*10* 

8.9785*10* 

9.4900*16* 

1.0010*11  * 

1  .0540*11* 

1.1379*11* 

1  .1630*11* 

1  .2191 *1 1 . 

1.2764*11* 

1.3351*11* 

1.3351*11  * 

1.4566*11  * 

1.5198*11, 

1.5846*11* 

1.6514*11* 

1.7202*11* 

1.7912*11  * 

1.3647*11, 

1.9408*11* 

2.0201* 11* 

2.1626*11 » 

2.1690*11* 

2.2796*11. 

2.3756*11* 

2.9781*11* 

3.3379*11 

2.4775*11* 

3.1472*11* 

2.5862*11* 

3.3615*11* 

2.7339*11* 

2.8333*1 1. 

▼<101  - 

3.0785*05* 

7.3785*35* 

7.2362*05* 

7.1927*35* 

7  .3979*35. 

7.0014*05* 

6. **037*05. 

6.8646*05* 

6.7039*05* 

6.601 7*05. 

6.4975*05* 

6.3917*35* 

6.2840*05* 

6.1  74  5  *35* 

5  .3631  *35. 

5.  94  96*05* 

5.8341*05* 

5.7164*05* 

5.5965*05, 

5.4743*05. 

3.3497*05* 

5.2228*05* 

5.0330*05* 

4.9525*35, 

4  .8254*35, 
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TABLE  10,  EOSTAB  INPUT 


BERYLLIUM  (Continued) 


9.; 9?5*os. 

9.5963*35* 

9,9319*35. 

9.2595*35, 

9.1033*05, 

3. 9983*05, 

3.7892*05, 

3.6  263*05. 

3.9538*05, 

3.2865*05, 

3.1 J33*05. 
2.1868*05, 
1.7191*05 

2.9269*05. 

1.9966*05* 

2.7191*05. 

1.7619*05, 

2.5933*35, 

2.3983*35, 

V501  = 

3.3323  « 

1.7193-  37* 

8.6269-07, 

3. 261.-36. 

9.9239-36, 

2. 5552-05. 

5.7596-05* 

1.1671-09. 

2.1685-09, 

3. 7511-C9, 

3.1127-09. 

9.9753-39, 

1.9378-33. 

2.0175-33, 

2  .8030-03, 

3. 7915-03. 

5 . Cl 10-03* 

6.9901-03, 

8.2582-03. 

1.0396-02, 

1.2795-02. 

1.5608-32. 

1  .839  9-  32. 

2.L J73-32, 

2.6795-02, 

3. 1989-02. 

3.6823-02. 

9.2807-02, 

9.9E09-C2. 

5.6991-02, 

6.5359-02. 

7  .9539*  02  . 

3.51 20-02. 

9.67  37-'2. 

1.0990 -Cl, 

1. 2969-01. 
2.9078-01 . 
9. "261-01 

1. 9196-01. 
2.8113-31  , 

1.6C69-C1, 
3.37  31-01  * 

1 .8227-01. 

2.09C6-01 , 

f  631  = 

3. 3591.11. 

3.8779*11  . 

3.7178* 11. 

3.7592*11. 

3.7985*11, 

3.8339.11  * 

3.3732*11 , 

3.91  99*11. 

3. "5*9*11* 

J .99 "2  *1 1  , 

9.0368.11. 

9.0780* 11. 

9.1167*11, 

9.1550*11. 

9.1925*11, 

9.2299*11. 

9.2659*11, 

9.3335*11, 

9.33*6*11, 

8.3676.11, 

9.3993*11. 

9.9296*11, 

9.9585*11, 

9.9857*11, 

9.5112*11, 

9.5399*11* 

9.5565*11. 

9.57  53*11* 

9.5927  *1  1  . 

9  .6070*1  1  * 

9.6189*11* 

8 • 8  8  3" *1 1 » 
3.8379*11 

8.6265*11. 

9.3596*11. 

9.6310*11. 

9.2?CO«ll. 

9.6315*11* 

9.(279*11, 

V70|  r 

5.3710*09, 

1.7529*35. 

1.83  79*05. 

1 .9298*35. 

2.3372*05. 

2.0866*05, 

2.1628*05, 

2.2259*05, 

2.3061*05, 

2.3732*05, 

2.9379*05* 

2.9935*35* 

2.55  >9*05, 

2.6111*05. 

2.5525*35, 

?• 7105*05. 

2.7550*05, 

2.7959*05, 

2.8332*05, 

2.8668*05, 

' 

2.8966*05* 

2.9225*05, 

2.91 85*05, 

2.9629*35. 

2  .9762*05, 

2.9858*05, 

2. "909*05. 

2.9915*05, 

2.9873*05, 

2.9  78?  *05. 

2.9637*05* 

2.9938*35. 

2.91 77*05, 

2.8352*35, 

2.8956*05, 

2. 7980*05. 
2.3832*05* 
1. 7191*05 

2.7915*05, 

2.2931*35, 

2.6  79  7*05, 
2.3599*"5» 

2.5955*05, 

2.5019*05, 

7  831  = 

9.6156*06* 

1.1997*07, 

-3.5629*11,- 

3.6  706*T6 » 

1.9282*11, 

-1.7933*11*- 

3.8695*11* 

9.93  97*11,- 

1.3538*11, 

1.9939*10 

21 


TABLE  11 


EOSTAB  INPUT :  TITANIUM 

M  = 

3.0000*02* 

3.2400*03* 

3.4  £59*03* 

3 • 73C0*C3 

4.2230*03* 

4.4653*03* 

4.7103*03* 

4.9553*33 

5.4450*03* 

5.6900*03* 

5.9350*03* 

6.1800*03 

6.5 730*03* 

6.9153*33* 

7.1*00*03* 

7.4353*33 

7*  8950*03 • 

8.1400*03* 

8.3850*03* 

8.6300*03 

9.1230*03* 

9.3650*33* 

9.6133*03* 

9.8550*33 

1.0745  *04  • 

1.0540*04* 

1.0835*04* 

1.1C8C*C4 

1.1570*04* 

1.2550*04 

1.1815*04* 

1.2360*04* 

1  .2315*34 

mi; 

0.0000  • 

2.8880*03* 

1.1178*04* 

3.6741*04 

2.5396*05* 

6.0264*35* 

1.2645*06* 

2.4636*36 

7.8730*06* 

1.3061*07* 

2.0795*07* 

3.1934*07 

5.8681*07* 

9.6826*37* 

1.3146*08* 

1.8025*28 

3.1176*08* 

4.0065*03* 

5.0785*08* 

6.3590*08 

9.5461*08* 

1  .1727*39* 

1 *4J  37*09* 

1.6814*39 

2.  3438*0  9* 

2.7404*09* 

3.1854*09* 

3.6823*09 

4  .8474*09* 
7.4800*09 

5.5233*39* 

6.2*83*09* 

7.0856*39 

f  2  3 1  = 

4.5100*00* 

4.1080*00* 

4.0703*00* 

4.031  8*PC 

3.9525*00* 

3.9115*30* 

3.8! 97*03* 

3.8259*33 

3. 7382*00* 

3.6923*00* 

3.6452*00* 

3.5969*00 

3.4964*00* 

3.4441*33* 

3.3433*03* 

3.3353*30 

3.2190*00* 

3.1584*00* 

3.0956*00* 

3.0307*00 

2.8936*00* 

2.8213*03* 

2.7434*33* 

2 .6  6e4 *3  3 

2.4966*00* 

2.4047*00* 

2.3C7l*00* 

2.2028*00 

1.9666*03* 

1.0163*00 

1.8285*03* 

1.6*79*00* 

1  .4652*13 

f  331- 

0.0000  * 

1.6421*10* 

1.7891*10* 

1,9381*10 

2.2425*10* 

2.3932*13* 

2.5364*13* 

2.7172*10 

3.04  76*10* 

3.2176*10* 

3.3909*10. 

3.5680*10 

3.9342*13* 

4.1243*13* 

4.31 36*13* 

4.5184*13 

4.9357*10* 

5.154  0*  10* 

5.3797*10* 

5. R1 32*10 

6.1071*13* 

6.3695*13* 

6.6l35*lC* 

5.9338*13 

7.5522*10* 

7.4911*10* 

8.2530*10* 

8.6426*13 

T43i  = 

4.5353*13* 
1. 3224*1 1 

1.0064*11* 

1.3*84*11* 

1.1473*11 

4.6947*05* 

4.0502*05* 

3.9924*05* 

3.9338*05 

3.8143*05* 

3.7532*35* 

3.6313*35* 

3.6233*35 

3.4994*05* 

3.4  332*05* 

3.3(59*05* 

3.2973*05 

3.1552*05* 

3.0335*35* 

3.0393*35* 

2.9336*35 

2.7767*05* 

2.6957*05* 

2.6  125*05  . 

2.5272*05 

2.3498*05* 

2.2573*05* 

2  .1  621*35* 

2  .0633*35 

1.8572*05* 

1.7482*05* 

1.6  750*05 * 

1.51 71*05 

1 .2645*05* 

1.1239*25* 

9.3725*04* 

8.4343*34 

3.975f'»P3» 
5.200303* 
6.4250*P3t 
7 .853303* 
8.8750*03* 
1 .3103*04* 
1.1325*04* 


1.0452*05* 
4.5287*06. 
4.7501*07* 
2.3931*08* 
7.8723*06* 
1 .991009» 
4.2350*09* 


3.992C*0C* 

3.7331*00* 

3.5473*Pr* 

3.2779*30* 

2.9634*00* 

2.5836*00. 

2.0901*00* 


2.0892*10*  , 
2.8909*1 0. 
3.749f'*:n, 

4 .7240*1 C» 
5.8554*10* 
7.2332*10* 
9.0670*10* 


3.8744*05* 
3.5544*25* 
3.2274*05* 
2  .3553  *05* 
2.4397*05* 
1  .9623*05* 
1.3937*05* 
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TABLE  11,  EOSTAB  IN1UT :  TITANIUM  (Continued) 


7.507  5O4 

Y5C1  = 

3.3333  » 

4  . 6331-  37. 

1.7332-06* 

5.2471-35* 

1  .4313-05* 

3. 3316-05 » 

7.1928-05* 

1.4  313-04  * 

2.65  71-04, 

4.6483-04 , 

7.7254-04* 

1.2282-33* 

1.8734-03* 

2.7773-33* 

3.9862-03* 

5. 5739-C 3* 

7.6162-03* 

1.0197-02* 

1.34C7-C2* 

1.7344-02* 

2.2121-32  . 

2.7862-32* 

3.4390*32* 

4 .2737-32* 

5.2257-02* 

6. 3361-02* 

7.6294-02* 

9.1338-02* 

1.0880-01* 

1.290  5- Cl » 

1.52:0-01  * 

1.3034-31* 

2.1319-31* 

2.5313-31. 

2  .9541  -01  • 

3.5011-01* 

1.3163*00 

4.1909-01* 

5.0911-01* 

6.4012-01, 

Y6"lr 

1.1310*11* 

1.2847*11* 

1.2975*11* 

1 .3131 *11. 

1  .3232*1  1  , 

1.3360*1 1* 

1.3488*11* 

1.3615*11* 

1.3743*11* 

1.3870*11* 

1.3997*11  * 

1.4123*11  • 

1.4247*11* 

1.4371  *1  1. 

1  .4493*1  1  • 

1. 4611*1 1* 

1.4731*11* 

1.4 647*11 . 

1 .4  9r  b* 11 , 

1.5066*11* 

1.5172*11  * 

1.5272*11* 

1  .5366*11* 

1.5455*11* 

1.5535*1 1* 

1.5610*11* 

1.5674*11* 

1.5729*11* 

1.5772*11. 

1.5802*11* 

1.5817*11* 

1.5813*11* 

1.5739*11* 

1  .5737  *1  1  * 

1.5654*11* 

1.5528*11* 

1.3225*11 

1.5343*11* 

1.5  C6  9* 11  * 

1.4628*11* 

Y7C1- 

2.3730*04* 

9.1376*34* 

9.4236*04* 

9.7543*34* 

1.3359*35* 

1.0374*05* 

1.0669*05* 

1.0955*05* 

1 .1252*05* 

1.14  93*05, 

1.1757*05* 

1.-JJ405* 

1  .2740*35* 

1 .2465*35* 

1.2675*05, 

1.2874*05* 

1.3057*05* 

1.3224*05* 

1.3375*05* 

1.3508*05* 

1.3621*05* 

1.3715*35* 

1.3737*05* 

1.38J7*;5* 

1  .366?*05» 

1.3862*05* 

1.3834*05* 

1.3  775*35. 

1.3684*05* 

1.3558*03* 

1 .3392*35* 

1.3131 *05* 

1 .2919*05. 

1.2599*35, 

1  .2207  *"5* 

1.1729*05* 

7.5075*04 

1.1136*05* 

1.0286*05* 

9.3776*04* 

YBClr 

1.8770*06* 

2.1174*35.- 

4.97  54*10.- 

2.6615*35* 

1  .0943*1  0. 

5.23CC*06* 

1.1153* 11* 

4 . 0  CCO*  C 1 . 

o.ooco  , 

o.cccc  * 

•  4.0327*13*- 

7.7051*13* 

3.4793*13*- 

4.7472*39* 

2  .7844*08 
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GRAY  EQUATION  OF  STATE  PARAMETERS  FOR  ALUMINUM 
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GRAY  EQUATION  OF  STATE  PARAMETERS  FOR  BERYLLIUM 
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3.  RESULTS  AND  DISCUSSION 


The  eighteen  RIP  calculations  performed  for  this 
research  are  assigned  calculation  numbers  and  described  in 
terms  of  the  equation  of  state  model,  the  material,  and  the 
spectrum  in  Table  15  .  The  calculated  impulses  and  peak 
stresses  at  the  end  of  the  source  time  for  aluminum,  beryllium, 
and  titanium  are  presented  in  Tables  16,  17,  and  18  respectively. 
Note  that  in  most  cases  the  agreement  between  the  three  equa¬ 
tion  of  state  models  is  fairly  good.  It  should  be  noted  that 
the  calculated  data  reveals  no  obvious  trends  concerning  the 
tendency  for  one  equation  of  state  to  consistently  predict 
stresses  and  impulses  which  lie  above  or  below  the  predictions 
of  the  other  two  models. 
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TABLE  15 


SUMMARY  OF  RIP  CALCULATIONS 


Calculation  No. 

Equation  of  State 

Material 

Spectrum 

1 

EOS1 

Aluminum 

SPEC  1 

2 

EOS1 

Aluminum 

SPEC  2 

3 

EOS1 

Beryllium 

SPEC  1 

4 

EOS1 

Beryllium 

SPEC  2 

5 

EOS1 

Titanium 

SPEC  1 

6 

EOS1 

Titanium 

SPEC  2 

7 

EOS  2 

Aluminum 

SPEC  1 

8 

EOS  2 

Aluminum 

SPEC  2 

9 

EOS2 

Beryllium 

SPEC  1 

10 

EOS  2 

Beryllium 

SPEC  2 

11 

EOS  2 

Titanium 

SPEC  1 

12 

EOS  2 

Titanium 

SPEC  2 

13 

EOS  3 

Aluminum 

SPEC  1 

14 

EOS  3 

Aluminum 

SPEC  2 

15 

EOS  3 

Beryllium 

SPEC  1 

16 

EOS  3 

Beryllium 

SPEC  2 

17 

EOS  3 

Titanium 

SPEC  1 

18 

EOS  3 

Titanium 

SPEC  2 
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TABLE  17 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  modified  PUFF  equation  of  state  is  easy  to  use  and 
gives  results  which  are  comparable  to  the  more  sophisticated 
models;  this  feature  of  course  is  the  primary  reason  that 
modified  PUFF  is  so  widely  applied  in  hydrodynamic  calculations. 
However,  this  model  has  several  disadvantages  which  are 
inherent  due  to  its  simplicity: 

1.  It  is  highly  empirical  and  is  therefore 
dependent  on  UGT  or  electron  beam  gen¬ 
erated  data. 

2.  The  melting  and  vaporization  phase  tran¬ 
sitions  are  described  by  energy  levels 
which  have  no  volume  (or  pressure)  depen¬ 
dence  . 

3.  The  model  contains  no  explicit  caloric 
equation  of  state;  the  temperature,  which 
may  be  required  for  some  applications,  is 
therefore  unknown. 


The  P.IP  mixed  phase  equation  of  state  was  developed  in 
order  to  partially  relax  the  above  restrictions.  However,  we 
find  that  additional  work  is  required  before  that  model  can 
be  considered  a  useful  calculational  tool.  The  following 
improvements  are  recommended: 

1.  The  lowest  temperature-volume  point  on  the 
mixed  phase  boundary  (see  Figure  3)  should 
be  the  boiling  point  at  ambient  pressure, 
i.e.,  0  =  0]-,,  the  boiling  temperature  and 
V  =  Vj3 ,  the  specific  volume  at  the  boiling 
temperature . 

2.  The  join  function  should  extend  from  the 
ambient  volume,  VQ,  to  the  boiling  volume, 

Vi.  (see  Figure  3)  .  Thus  the  join  function 
will  not  penetrate  the  mixed  phase  region 
causing  the  discontinuities  which  exist  in 
the  present  RIP  mixed  phase  equation  of 
state . 
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Pressure 


Join  function 
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3.  The  EOSGEN  code  should  be  incorporated  into 
the  RIP  code  as  a  subroutine  which  is  called 
after  the  initial  call  to  the  EOS2  subroutine; 
this  would  eliminate  the  need  to  perform  a 
separate  calculation  with  the  EOSGEN  code  and 
would  greatly  simplify  data  input  for  the  RIP 
mixed  phase  equation  of  state. 

The  GRAY  equation  of  state  contains  the  most  sophisti¬ 
cated  equation  of  state  models  available  in  the  RIP  code.  How¬ 
ever,  the  procedure  used  to  join  the  Grover  liquid  equation  of 
state  and  the  Young-Alder  model  destroys  the  physical  signi¬ 
ficance  of  the  Young-Alder  model,  particularly  in  the  liquid- 
vapor  region.  Additionally,  since  the  join  volume  used  in  the 
GRAY  equation  of  state  is  too  large,  the  Grover  model  is 
applied  in  regions  where  its  accuracy  is  highly  questionable. 

The  procedure  we  recommend  for  correcting  these  deficien¬ 
cies  is  best  explained  in  terms  of  a  pressure-volume  diagram 
which  is  presented  schematically  in  Figure  4.  The  reference  or 
ambient  condition  is  assumed  to  be 

P  =  Pq  =  0  ,  pressure 

T  =  Tq  ,  temperature 

and  V  =  VQ  ,  specific  volume 

The  specific  volumes  indicated  on  the  abscissa  are: 

VmQ  =  reference  pressure  intercept  of  solidus 

VLq  =  reference  pressure  intercept  of  liquidus 

V.  =  reference  pressure  intercept  of  liquid- 
vapor  phase  boundary. 

V  =  critical  specific  volume 
c 
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■olidu*  y  liquidu* 


Isotherms  passing  through  the  reference  pressure  at  Vq,  Vmo» 
Vh0  and  Vc  are  denoted  by  To,  Tmo/  Tbo  and  Tc  respectively. 

The  basic  problem  with  the  GRAY  equation  of  state  is  that  the 
join  volume  lies  to  the  right  of  Vbo  in  Figure  4;  thus  the 
Young-Alder  isotherms  in  the  liquid-vapor  region  are  badly 
distorted.  Consequently,  no  attempt  was  made  to  perform  a 
Maxwell  construction  in  order  to  identify  the  liquid-vapor 
phase  boundary;  failure  to  perform  this  construction  results 
in  the  appearance  of  unstable  states  to  the  right  of  the  join 
volume  due  to  loops  in  the  Gibbs  free  energy  along  isotherms. 
We  therefore  recommend  that  the  Grover  and  Young-Alder  models 
be  joined  at  some  locus  between  the  liquidus  and  Vfc0,  the 
specific  volume  at  the  reference  pressure  and  the  boiling  tem¬ 
perature.  Joining  the  two  models  in  this  region  will  leave 
the  Young-Alder  model  in  tact  and  allow  one  to  perform  the 
Maxwell  construction.  One  of  several  methods  may  be  used  to 
join  the  Grover  and  Young  Alder  models: 

1.  The  models  could  be  joined  at  the  locus  of 
intersection  of  corresponding  isotherms  if 
that  locus  does  not  lie  within  the  mixed 
phase  region. 

2.  The  two  models  could  be  joined  by  a  fairing 
function  applied  between  VLq  and  VbQ. 

3.  A  fairing  function  may  be  applied  to  the 
entire  region  between  the  liquidus  and 


It  should  be  emphasized  that  the  recommended  changes 
to  the  GRAY  equation  of  state  are  very  important  because  they 
will  make  it  possible  to  determine  all  phase  transitions 
explicitly;  this  capability  is  especially  desirable  for  appli 
cations  involving  the  x-ray  induced  response  of  porous  mate¬ 
rials  .  A  recent  study  performed  at  Lawrence  Livermore  Labor¬ 
atory  I  HI  shows  that  impulse  production  in  porous  materials 
occurs  on  two  distinct  time  scales.  The  short  time  contri¬ 
bution,  which  occurs  in  less  than  a  microsecond,  is  due  to 
pore  closure  by  thermal  expansion  which  results  in  high  pres¬ 
sures  and  substantial  blow-off  velocities.  The  second  con¬ 
tribution,  which  occurs  over  a  time  period  of  tens  of  micro¬ 
seconds,  is  due  to  the  vapor  pressure  in  the  mixed  liquid/ 
vapor  region.  The  vapor  pressure  tends  to  accelerate  the 
material  in  front  of  it  (either  blowoff  material  or  material 
in  the  missile  structure  which  acts  as  a  tamper)  over  a  long 
time  period  and  can  result  in  a  significant  contribution  to 
the  total  impulse.* 


It  should  be  mentioned  that  explicit  identification  of  phase 
boundaries  is  also  important  with  respect  to  laser  induced 
material  response  because  it  would  enable  one  to  identify 
regions  in  which  liquid  droplets  are  present.  Liquid  droplets 
can  drastically  affect  the  optical  properties  of  the  medium 
and  these  properties  are,  of  course,  quite  important  for  pre¬ 
dicting  laser  interactions . 
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APPENDIX  A 
THE  EOSGEN  CODE 


Three  tasks  are  performed  by  the  EOSGEN  code:  1)  Eval¬ 
uation  of  the  seven  coefficients  in  the  virial  expansion, 

2)  definition  of  the  mixed  phase  boundary,  3)  evaluation  of 
the  seven  coefficients  in  the  join  function.  These  tasks  are 
discussed  separately  in  the  following  sections  of  this  Appen¬ 
dix. 


A. 1  THE  VIRIAL  EXPANSION 


The  virial  expansion  employed  by  the  RIP  mixed  phase 
equation  of  state,  with  density  and  energy  as  independent 
variables,  is 


P  = 


(zlP+z2p2+z4p3jE 


zllp+z12p  +z13p 


+  z14p4+z15p5 


(1) 


where 


3P  V  9 


Zi  = 


3P  V 
c  c 


i  e 


Z2  = 


-3P  V“+V  [b  V  +  2R  ( 0  )1 
c  c  oL  o  o  oj 

9  -0 
c  o 


-V  0  |e  V  +2R { 0  )1  +  3P  V20 
_  o  clo  o  o  J _ c  c  o 


9-0 
c  o 


Preceding  page  blank 
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In  Equation  (1),  P  is  the  pressure,  p  is  the  density  and 
E  is  the  energy.  The  coefficients,  z±,  are  expressed  in  terms 
of  the  critical  and  reference  point  characteristics  presented 
in  Tables  4,  5  and  6.  The  coefficients  z\  through  Z7  are 
evaluated  by  the  EOSGEN  code  and  must  be  read  into  the  first 


i 

j 
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seven  locations  of  the  Z  array  in  EOSTAB  common  in  the  RIP 
code.  The  coefficients  zn  through  zi5  are  evaluated  by  sub¬ 
routine  E0S2  in  the  RIP  code  and  need  not  be  input. 


A .  2  DEFINITION  OF  THE  MIXED  PHASE  BOUNDARY 


Corresponding  values  of  pressure,  density,  energy  and 
sound  speed  on  the  mixed  phase  boundary  are  required  as  input 
for  the  RIP  mixed  phase  equation  of  state.  These  data  are 
determined  by  locating  points  having  equal  pressure  and  Gibbs 
free  energy  on  a  series  of  isotherms  between  the  reference  tem¬ 
perature,  80,  and  the  critical  temperature  6C.  Referring  to 
the  schematic  pressure/volume  diagram  presented  in  Figure  A.l, 
the  two  points  on  the  mixed  phase  boundary  of  an  isotherm  with 
temperature,  0O  <  0  <  0C,  are  determined  by  the  following 
equalities: 

Pu(6'  pu}  =  PL<0'  pL>  *  (2) 

Gu(0,  Pu)  =  Gl(0,  Pl)  ,  (3) 

where  G  is  the  Gibbs  free  energy,  P  is  the  pressure  and  pl  and 
pu  are  densities  on  the  low  and  high  density  sides  of  the  mixed 
phase  boundary  respectively.  The  pressure,  as  a  function  of 
density  and  temperature,  is  given  by 

P  (p ,  8)  =  Zj0p  +  (z20+z3)p2  +  (z40+z5)p3  ,  (4) 

and  the  Gibbs  free  energy  is  defined  as 
G  =  E+PV-0S  , 


where  E  is  the  specific  energy  and  S  is  entropy.  Taking  the 
differential  of  the  Gibbs  free  energy  yields 

dG  =  dE+PdV+VdP  -  0dS-Sd0  , 

which  reduces  to 

dG  =  VdP  (5) 
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Pressure 


Figure  A.l  Intersection  of  isotherms  with  mixed  phase  boundary. 
(Note:  PL  =  1/VL;  Pu  =  1/VU)  . 
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for  an  isothermal  process.  Therefore,  the  Gibbs  free  energy, 
referenced  to  the  initial  density,  is 


Substituting  Equation  (4)  into  the  above  relation  and  perform¬ 
ing  the  indicated  integration  yields 

G(e,P)  =  -*1e  Jn^2)  +(26*.,  +»3)(p-p0)  +  f(z40«5)(p2-p2)(.^ 

Now,  Equations  (4)  and  (7)  are  substituted  into  Equation  (2) 
and  (3)  respectively,  and  the  quantities  Pu>  Pu  an d  pl  are 
calculated  by  a  numerical  iteration  scheme  in  the  EOSGEN  code. 
The  specific  energies  on  the  mixed  phase  boundary  are  deter¬ 
mined  by 


E 


u 


E(9,pu)  , 


and 


El  =  E(6,pl)  , 


where 


E 


z3p 


5  2 

+  P  +  Zg0+z^ 


I 


is  the  caloric  equation  of  state  employed  in  the  expanded 
region  of  the  RIP  mixed  phase  equation  of  state.  The  sound 
speeds  on  the  mixed  phase  boundary  are  given  by 


evaluated  at  0  and  P  =  Pu 
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evaluated  at  0  and  p  =  p 


c 


L 


L  ' 


The  quantities  on  the  mixed  phase  boundary  are  stored 
in  the  following  arrays  in  EOSTAB  common  in  the  RIP  code: 


0  THETA  (I) 
Pu  ■>  PU(I) 

Pu  +  RHOU(I) 
Eu  +  EU (I ) 
cu  -  CU(I) 

PL  +  RHOL(I) 
El  ->  EL  (I) 

CL  +  CL (I)  . 


The  number  of  data  points  to  be  stored  in  each  of  the 
above  arrays  must  be  input  to  Z  (8)  of  the  Z(I)  array  which  is 
also  in  EOSTAB  common. 


A. 3  THE  JOIN  FUNCTION 


The  purpose  of  the  join  function  is  to  allow  pressure 
and  sound  speed  continuity  between  the  Mie  Gruneisen  solid 
equation  of  state  for  the  compressed  region  and  the  virial 
expansion  for  the  expanded  region.  The  join  function  is  applied 
over  a  density  range  beginning  at  the  ambient  density,  p0, 
and  ending  at  a  slightly  expanded  cut-off  density  of  pco 
(see  Figure  A. 2).  It  was  determined  that  in  general  it  is 
possible  for  the  join  function  originally  reported  in  Refer¬ 
ence  1  to  contain  several  singularities  over  this  density 
range;  in  fact,  these  singularities  appeared  when  the  calcula¬ 
tions  performed  for  this  report  were  first  attempted.  We 
therefore  decided  to  develop  a  new  join  function  which  would 
not  exhibit  singularities  under  any  conditions. 
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Expanded  Region 
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Figure  A. 2  Schematic  representation  of  the  RIP  mixed  phase  equation  of  state. 


The  new  join  function  developed  for  this  series  of  cal¬ 
culations  is : 


c  +c  1  J^-j 

5  6\po-poo  /  7VVpoc / 


C1+C2P+C3P  +C4P  ' 


(7) 


where  P2  is  the  pressure  and  the  subscript,  2,  serves  to  iden¬ 
tify  the  join  function.  The  seven  coefficients,  ci***C7,  are 
determined  by  continuity  conditions  at  the  ambient  and  cut-off 
densities.  Letting  the  subscripts  1  and  3  represent  the 
Mie  Gruneisen  solid  equation  of  state  and  the  virial  expansion 
respectively,  we  have 


P1  =  GoPoE  +  (  A1M+A2y2+A3y3)  (l  -  , 


(8) 


where 


y  = 


t'1  • 


P-,  = 


ZXP+Z2P  +Z4P' 


+  z11P+z12p  +zi3p  +Z14P  +zi5p" 


(9) 


The  appropriate  boundary  conditions  for  determining  the  coeffi¬ 
cients  in  the  join  function  are: 


pl(po,E)  -  p2(po-E) 


(10) 


3P. 


3P, 

w 


at  p  =  p. 


(11) 
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9E 


9P, 


9E 


at  p  =  pQ  , 


P2(p=o'E)  -  P3(»co'E) 


(12) 

(13) 


3P2 

TT 


9P. 


E 


9  P 


at  P  =  pCO  • 


(14) 


9P. 

9F 


9P; 

9E" 


at  p  =  p 


co 


(15) 


Substituting  Equations  (7) ,  (8) ,  and  (9)  into  boundary  condi¬ 
tions  (10)  through  (15)  yields  the  following  set  of  indepen¬ 
dent  equations  for  the  undetermined  coefficients  in  the  join 
function. 


c5  =  N1 


(16) 


C5+c6+c7  =  N4 


(17) 


2  _  .  3 

-  c,  +  -  c_  =  -N  , 

P^-P^  6  P^-P.  -  7  6 


O  CO 


O  OC 


(18) 


2  3 

c.+c0p  +c0p  +c.p  =  0 
1  2o  3  o  4  o 


c2po+2c3po+3c4po  =  N2 
Cl+C2pco+C3pco+C4pco  =  N3 


V2c3pco+3c4pco  =  N‘ 


(19) 

(20) 
(21) 

(22) 


where 


N1  =  poGo 


N2  =  Ax 
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N.,  =  Z,  ,  p  +Z,  ~p  +Z,,P  +Z,  *  P  +Z,rP 

3  llHco  12Hco  13Hco  14Hco  15Hco 


N4  = 


2  3 

Z,  p  +Z-.P  +z.p 
lMco  2hco  4kco 


N5  =  zU+2z12pco+3z13pco+4z14pco+5z15pM 


N6  - 


z.+2z^p„  +3z . p 
1  2Kco  4  co 


Equations  (16)  through  (21)  are  easily  solved  by  Cramer's  rule; 
the  task  is  performed  by  the  EOSGEN  code.  Additionally,  the 
cut-off  density,  pco,  is  selected  automatically  by  the  EOSGEN 
code  based  on  the  following  considerations : 

1.  The  coefficient  of  E  in  Equation  (9)  has 
a  value  at  p  =  pCq  which  is  less  than  the 
coefficient  of  E  in  Equation  (8) . 

2.  The  remaining  terms  in  Equation  (9)  have  a 
negative  value  at  p  =  pco  and  their  deriva¬ 
tive  has  a  positive  value  at  p  =  p  . 

The  new  join  function  discussed  in  the  above  paragraphs 
is  mathematically  cleaner  (no  singularities)  and  greatly  facil¬ 
itates  application  of  the  RIP  Mixed  Phase  equation  of  state. 


48 


APPENDIX  B 
CORRECTIONS  TO  EOS 2 


i 


Two  types  of  corrections  to  subroutine  EOS2  (the  RIP 
mixed  phase  equation  of  state)  were  required  in  order  to  per¬ 
form  the  calculations  presented  in  this  report.  The  first 
set  of  corrections  are  related  to  the  incorporation  of  the 
new  join  function  (see  Appendix  A. 3)  into  subroutine  E0S2. 

The  second  correction  is  a  temporary  measure  designed  to  com¬ 
pensate  for  inherent  discontinuities  in  the  RIP  mixed  phase 
equation  of  state.  These  discontinuities  exist  between  the 
mixed  phase  region  and  the  virial  expansion  region  below  the 
zero  pressure  line  in  Figure  B.l;  the  mixed  phase  and  virial 
expansion  regions  in  question  are  denoted  as  27  and  26, 
respectively,  in  the  same  figure.  The  task  of  correcting  the 
discontinuity  in  a  thermodynamically  consistent  manner  is 
beyond  the  scope  and  intent  of  this  program;  however,  a  tem¬ 
porary  correction  was  required  in  order  to  avoid  numerical 
convergence  problems  encountered  when  a  Lagrangian  cell 
attempted  to  follow  a  path  from  the  26  region  to  the  27  region. 
Since  such  a  transition  is  discontinuous  and  physically  meaning¬ 
less  and  because  convergence  problems  were  encountered  whenever 
such  a  transition  occurred,  it  was  decided  to  artificially  for¬ 
bid  it  with  appropriate  modifications  to  the  E0S2  subroutine. 

A  FORTRAN  listing  which  details  the  required  correc¬ 
tions  to  the  EOS 2  subroutine  in  the  RIP  code  is  presented  at 
the  end  of  this  appendix. 
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The  FOLLOWING  CHANGES  to  EOS2  I r I p  mixed  PHASE  equation  OF  STATE)  ARE 

Required  in  the  existing  model  in  order  to  incorporate  the  join  function 

AND  TO  COMPENSATE  FOR  D I  SCOUT  I NU I TES  BETWEEN  THE  MIXED  PHASE  REGION  AND 
The  VIKIAL  EXPANSION  region  below  the  zero  pressure  line. 


1.  IN  THE  CODING  WHICH  FOLLOWS  •••  POINT  IS  IN  MIXED  PHASE  REGION  ••• 
COMMENT  CARD  ADD  The  FOLLOWING  LINES  AFTER  THE  LINE  WHICH  STATES! 

TEMPI  1  I ) »PU ( NM A  X ) 

add  lines: 

IF(TEMPII)  .GT.  TEMPI  1 C )  )  GO  TO  155 
P I J I =  TEMP ( 10) 

I  VFLaGI J>»26 
Cl J)«ESTC0N(6,N) 

GO  TO  260 

AND  CHaNGE  THE  STATEMENT  WHICH  STATES  t 

I  F  (  TEMP (  1  )  .GT.  TEMP  ill))  TEMPI  1 ) “TEMP  I  1  I  ) 

to  the  following: 

1SS  IFITeMPII)  .GT.  TEMP  I  I  1  )  )  TEMP  I  1 ) -TEMP  I  l  1  ) 


2.  AFTER  THE  LINE  OF  CODING. 

?SC  I VFL AG  I J ) *25 

The  following  LINES  have  been  deleted  and  replaced: 

TEMP(2)»C0SIESTC0NI22,N)*TEMP{ | ) ) 

TEMP (3)«ESTcON  121 |N) / TEMP  I  2 > ♦ESTCON  I  23 , N ) 

TEMPIH)»S1NITEMPI 1 ) ) 

TEMP!S)»ESTC0N121  ,N)*ESTcON122,N)*TEMPi1«)/I  l  ESTCON  I  1  ,N) -ESTCON  I  16, 
,  N ) ) *TEMP I  2 ) *  *2 ) 

the  new  lines  should  read: 

TEMPI2)»TEMP« 1 ) • *2 

TEMPI3)»ESTC0N(23iN)*TEMPI2)*IESTC0N(21 , N ) ♦ESTCON I  22 , N ) ♦TEMP  I  1  ) ) 
TEMP  I 5)«- I  2. *EST CON  I  2 1 , N ) »TEMP 1 1 ) +3 . *ESTCON ( 22 , N ) • TEMP  I  2 ) ) / 

.  (ESTCONI 1 ,N)-ESTCON« 16.N) ) 


51 


APPENDIX  C 
CORRECTIONS  TO  EOS3 


Two  coding  errors  were  discovered  in  the  expression  for 
temperature  in  the  liquid-vapor  region  of  the  GRAY  equation  of 
state.  The  FORTRAN  changes  required  to  correct  these  errors 
are  presented  in  the  listing  at  the  end  of  this  appendix. 

Due  to  the  lack  of  continuity  between  the  liquid  vapor 
region  and  the  liquid,  melt,  and  solid  regions  (see  Figure  C.l) 
problems  (complex  temperatures)  were  encountered  whenever 
Lagrangian  cells  attempted  to  cross  the  boundary  between  them; 
this  problem  is  inherent  to  the  GRAY  equation  of  state  because 
the  liquid-vapor  region  is  continuous  only  with  the  hot  liquid 
region.  The  task  of  providing  complete  continuity  in  this 
equation  of  state  is  beyond  the  scope  and  intent  of  this 
project.  Therefore,  temporary  FORTRAN  changes  were  incorpor¬ 
ated  into  the  E0S3  and  EOS  subroutines  which  prevent  the 
liquid-vapor  region  from  being  entered  by  any  path  not  leav¬ 
ing  the  hot  liquid  region.  These  changes  are  detailed  in  the 
FORTRAN  listing  which  follows. 


Preceding  page  blank 
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Energy 


Figure  C.l  Schematic  representation  of  GRAY  EOS. 


54 


f  HE  FOLLOWING  CHANGES  TO  EOS  AND  EOS3  (GRAY  EQUATION  OF  STATE  >  SUBROUTINES 
ARE  RtQUIREO  IN  THE  EXISTING  MODELS. 

A.  SUBROUTINE  EOS 

I.  IN  THE  SECTION  OF  CODING  FOLLOWING  THE  COMMENT  STATEMENT  •••  LRL  EOS  ••• 
ADD  TWO  LINES  AFTER  ThE  ONE  WHICH  STATES  I 

30  RSAVE»RHO(J) 

ADDED  LINES! 

PSAVE2-PI J) 

ISAVE-IVFLAGI J) 

ALSO  IN  THIS  SAVE  SECTION  TWO  LINES  ARE  ADDED  IMMEDIaTLY  FOLLOWING  THE 

LINE  stating: 

P  S  A  V  E  ■«  P  (  J) 

ADDED  LINES: 

I VFLAGI J)*ISAVE 
P ( J ) -PSAVE2 


b.  SUBROUTINE  E0S3 

1.  IN  ThE  SECTION  OF  CODING  FOLLOWING  THE  COMMENT  STATEMENT  •••  VAPOR  REGION 
•••  add  Six  lines  OF  CODING  between  The  COMMENT  STATEMENT  and  The  line 
which  reads: 

90  IVFLaG<J)»3S 

90  IF(IVFLAG(JI  .GE.  33)  GO  TO  8S 
P( J)«0. 

GO  TO  too 

B5  I F ( P ( J  )  .GT.  I.E-12)  GO  TO  86 
P  (  J  )  *0  « 

GO  TO  100 

ANO  change  the  line  of  coding  which  reads: 

90  I VFL AG ( J  )  *35 

to  read: 

86  IVFLaG(J)>3S 


55 


ALSO  •  IN  THE  SAME  SECTION  OF  CODING  CHANGE  THE  STATEMENT  WHICH  ST  a  T  E  b 


TEMPI7  »«-IE(  Jl^ECONI  I7»N)*TEMP(  I  )/ESTC0N<26,N)-EC0N(*I,NI  *TEMPI&|- 
.  ECON(lOiN) 

TO  THE  FOLLOWING: 

TemP(7)»-(EIJ)>ECON(I7,N)*TEMP(1>/ESTCONI26,N) 

1  -EC  ON  I  |H,NI*TEMP(M-ECONI  I0,N)»2> 

ANu  ADD  the  following  is  lines. 

T1  2  3*I«TEMPI5  >  ••2-2*  «TtMPI  HI  •  TEMP  I  7  I 
IMTI2JM  ,GT.  J)  GO  TO  >2 

PRINT  *  I  .ECONI  S  ,N)  .ECONUUiNI  ,E  CONI  I  I  ,N  »  ,  EC  UN  |  I  2  ,  N  I  ,  t  C  ON  (  1  H  ,  N  )  , 

I  ECON  I  1 S ,N  I  ,ECOM I  6 ,N )  ,ECON( 1 7 ,N ) ,TEMP I  6  ) 

9 |  F0RM*TI5* .VhNPKIP1  *  ,  E  I  2  •  3  ,  /  /  , 

1  SX.SHD1  -  ,  E  I  2  •  J  ,  /  /  , 

2  S *  ,  S h D  2  >  t  E I  2  •  3 1  //  ■ 

3  S * . ShD 3  «  ,£12.3,//. 

H  SA.ShCI  ■  ,  E 12  »  3 , // , 

S  SA.ShC?  ■  ,  E 1 2  *  J , / / , 

4  SA.SHC3  *  ,£ I  2 . 3 , / / , 

7  SX ,  10HAYPKIME  *  , E 1 2 . 3  ,  /  /  , 

8  SA.SHFE  *  ,  E 1 2 • 3 , / / ) 

CALL  SPRINT  (AHEOS3  .AHLRR0R9) 

92  CONTINUE 
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